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Molecular absorption and photo-electron spectra can be efficiently predicted with real-time time-dependent 
density- functional theory (TDDFT). We show here how these techniques can be easily extended to study time- 
resolved pump-probe experiments in which a system response (absorption or electron emission) to a probe 
pulse, is measured in an excited state. This simulation tool helps to interpret the fast evolving attosecond 
time-resolved spectroscopic experiments, where the electronic motion must be followed at its natural time- 
scale. We show how the extra degrees of freedom (pump pulse duration, intensity, frequency, and time-delay), 
which are absent in a conventional steady state experiment, provide additional information about electronic 
structure and dynamics that improve a system characterization. As an extension of this approach, time- 
dependent 2D spectroscopies can also be simulated, in principle, for large-scale structures and extended 
systems. 



I. INTRODUCTION 

Pump-probe experiments are the preferred tech- 
nique to study the dynamical behaviour of atoms and 
molecules: the dynamics triggered by the pump pulse 
can be monitored by the time-dependent reaction of the 
system to the probe pulse, a reaction that can be mea- 
sured in terms of, for example, the absorption of the pulse 
intensity, or of the emission of electrons.^ The time reso- 
lution of these experiments is mainly limited by the du- 
ration of the pulses - although it is also limited by the 
ability of the experimenter to ascertain their relative time 
delay and shape. In order to precisely fix this delay, the 
two pulses are coherently synchronized - in fact, they 
have the same origin, or one of them is used to generate 
the other -, so that the delay is gauged by an optical path 
differenceP The electron dynamics has a natural time- 
scale in the range of attoseconds, and therefore could 
not be studied by pump-probe spectroscopy until the ad- 
vent of attosecond-pulse laser sources one decade agoP^ 
Nowadays time-resolved spectroscopy can be utilized to 
monitor electron dynamics in real time, giving birth to 
the field of attosecond physics.^ 

A wealth of possibilities exists, depending on the fre- 
quencies, durations and intensities of the two pulses. A 
common set-up in attosecond physics employs a XUV 
attosecond pulse and the relatively more intense NIR or 
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visible longer (few femtoseconds) pulse used for its gener- 
ation. Combining two XUV attosecond pulses is in prin- 
ciple possible (and has been theoretically analysed^) , but 
unfortunately the low outputs of current XUV attosec- 
ond pulses lead to much too weak signals. Another choice 
to make is the final observable, i.e. what kind of system 
reaction is to be measured as a function of the time de- 
lay. In this work we focus on two common choices. One, 
observing the emission of electrons (their energies, angu- 
lar distribution, or total yield) from the pumped system 
due to the probe pulse. This can be called time-resolved 
photo-electron spectroscopy (TRPES). Two, observing 
the optical absorption of the probe signal, which can be 
called time resolved absorption spectroscopy, or transient 
absorption spectroscopy (TAS). 

Both techniques can of course be used to look at longer 
time resolutions - and there is already a substantial body 
of literature describing such experiments. If we look at 
molecular reaction on the scale of tens or hundreds of 
femtoseconds, the atomic structure will have time to re- 
arrange. These techniques are thus mainly employed to 
observe modification, creation, or destruction of bonds, 
a field now named femtochemistry. 8 

TAS, for example, has been succesfully employed to 
watch the first photo-synthetic events in cholorophylls 
and carotenoids, 9 that transform the energy gained by 
light absorption into molecular rearrangements. A review 
describing the essentials of this technique can be found 
in Ref. [10J Note, however, that in addition to following 
chemical reactions, femtosecond- long pulses may also be 
used for example for characterizing the final electronic 
quantum state of ionized atomsP^ 

In TRPES, the probe pulse generates free electrons 
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through photo-ionization, and one measures their energy 
or angular distribution as a function of time; if this time is 
in the femtosecond scale one can follow molecular dynam- 
ics in the g as phase, as demonstrated already in the mid 
1990sp^ES although this technique had already been em- 
ployed to follow electronic dynamics on surfaces.^ This 
methodology is well documented in recent articles fHHUl 

If the goal -as in this work- is to study the electronic 
dynamics only, disentangling them from the vibronic de- 
grees of freedom, then one must move down these spec- 
troscopic methods to the attosecond regime!^ In this 
regime, both TAS and TRPES have recently been demon- 
strated. Regarding TAS, we may cite as a prototypical 
example the recent experiment of Holler et at^ where, 
the transient absorption of an attosecond pulse train (cre- 
ated by high harmonic generation) by a Helium gas tar- 
get, was studied in the presence of an intense IR pulse. 
The absorption was observed to oscillate as a function 
of the time-delay of pump and probe. Another example 
is the real-time observation of valence electron motion 
reported by Goulielmakis et al? 2 

Several cases of use of TRPES with attosecond pulses 
have also been recently reported. For example, Uiber- 
acker et a/P^ could observe in real time the light in- 
duced electron tunneling provoked by a strong NIR pulse, 
demonstrating how this electron tunneling can be used 
to probe short lived electronic states. Smirnova et a¥^ 
also studied the ionization of an atom by an attosecond 
XUV pulse in the presence of an intense laser pulse, as 
a function of the time delay between both. Jonnsson et 
aW^ employed attosecond pulse trains and a Helium tar- 
get, and attosecond photoelectron spectroscopy was also 
demonstrated to yield useful information for condensed 
matter systems.^ 

All these advances demand an appropriate theoretical 
modeling. The use of more than one pulse of light intrin- 
sically requires to go beyond any "linear spectroscopy" 
technique - although if the pulses are weak a perturbative 
treatment may still be in order. This non-linear behav- 
ior, provides much more information about the system 
at the cost of a increasingly difficult analysis. The use of 
two (or more) coherent pulses of light, with fine control 
over their shape (sometimes called a "multidimensional 
analysis"), permits a deeper characterization. This fact 
was already acknowledged in the field of nuclear magnetic 
resonance, or later in femtochemistry - see for example 
Refs. l27l and l28l for theoretical treatments of these cases. 

A recent theoretical analysis of attosecond TAS based 
on perturbation theory was given by Baggesen et a/P 
Gaarde et aP^ presented a study for relatively weak 
pumping IR pulses in combination with XUV ultrafast 
probes, for Helium targets and based on the single active 
electron approximation (SAE). Very recently, the exper- 
imented reported by Ott et al. 30 , in which the ultrafast 
TAS of Helium displayed characteristic beyond-SAE fea- 
tures, was theoretically analyzed in Ref. [3TJ utilizing an 
exact solution of the time-dependent Schrodinger equa- 
tion, that cannot however be easily extended to larger 



systems. Finally, the above-mentioned experiment of 
Goulielmalkis et aP^ was analized with the model de- 
scribed in Ref. [32j which treated the pump IR pulse non- 
perturbatively. 

Indeed, it would be desirable to analyze these processes 
with a non-perturbative theory (since at least one of the 
pulses is usually very intense), which at the same time 
is capable of going beyond the SAE and accounting for 
many-electron interaction effects. This last fact is rel- 
evant since the attosecond time resolution obtained in 
this type of experiments is able to unveil the fast dy- 
namical electron-electron interaction effects. The SAE, 
which essentially assumes that only one electron actively 
responds to the laser pulse, has been successfully used to 
interpret many strong-field processes. However, its range 
of validity is limited, and roughly speaking it is expected 
to fail whenever the energies of multielectron excitations 
become comparable to the laser frequencies or the single 
electron excitations) 3 - 3 -^ 

Time-dependent density functional theory (TDDFT^P^ 
in principle meets all requirements: it may be used non- 
perturbatively, includes the electron-electron interaction 
and can handle out-of-equilibrium situations. It has been 
routinely used in the past decades to study the electron 
dynamics in condensed matter in equilibrium. By this 
we mean that, usually, one computes the linear or non- 
linear response properties of systems in the ground state 
(or at thermal equlibrium). In pump-probe experiments, 
however, one must compute the response of a system that 
is being driven out of equilibrium by an initial pulse. In 
this work, we will explore the usability of TDDFT for 
this purpose, and show how, at least for the two cases of 
TAS and TRPES, the extension is straightforward. 



II. THEORY 

Density-functional theorjE' establishes a one-to-one 
correspondence between the ground-state density and 
the external potential of a many-electron system. This 
implies that any system property is, in principle, a 
ground-state density functional. The computation of 
the gro und-state density usually follows the Kohn-Sham 
(KS) 36 scheme, in which one utilizes a fictitious system 
of non-interacting electrons that has the same ground 
state density. For excited states, properties, however, 
or in order to simulate the behavior of the system in 
time-dependent external field s, one must use its time- 
dependent version, TDDFTl MMlM 

In the case of TDDFT, a one-to-one correspondence 
also exists between the time-dependent densities and 
potentials. One also uses an auxiliary fictitious sys- 
tem of non-interacting electrons that produces the same 
time-dependent density. This substitution is the source 
of the great computational simplification, since a non- 
interacting system of electrons can in general be repre- 
sented by a single Slater determinant formed by a set 
of "Kohn-Sham" orbitals, ^ (i = l,...,iV/2). We will 
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assume a spin-compensated system of N electrons dou- 
bly occupying N/2 spatial orbitals. If the real system 
is irradiated with an external field characterized by a 
scalar potential v(r,t) (the extension to vector poten- 
tials is also possible), the "time-dependent Kohn-Sham" 
(TDKS) equations that characterize the evolution of the 
fictitious system are (atomic units will be used hereafter): 

d 1 

i^i{r,t) = --V 2 <Pi(r,t) + v KS [n](r,t)(pi(r,t) , (1) 



N/2 

n(r,t)=2j2\Vi(r,t)\ 2 . 



(2) 



The time-dependent density n(f,i) is the central object, 
and is identical for the real and for the KS systems. The 
KS potential vks is a functional of this density, and is 
defined as: 

VKs[n](r,i) = v (r) +v(r,t) +m[n](r,t) +v xc [n](r,i) , 

(3) 

where the Hartree potential v^l corresponds to a classical 
electrostatic term 



vn[n](r,t) 



d 3 r / nftt) 



\f' — r\ ' 



(4) 



vo(r) is the static external potential that characterizes 
the system in its ground state (in a molecule, originated 
by a set of nuclei), and the "exchange and correlation" 
potential is v xc [n] . The exchange correlation potential is 
also a functional of the density and accounts for all the 
intricate many-electron effects. It is in practice unknown 
and must be approximated 0HEU 

The TDKS equations can be utilized, either directly or 
in appropriately transformed manner, to compute the re- 
sponse of a many-electron system to a perturbation, weak 
or strong. In the perturbative regime, ideally one wishes 
to obtain the response functions [(hyper)-polarizabilities, 
optical and magnetic susceptibilities, . . .], since (i) these 
objects then permit to predict any reaction in the ap- 
propriate order, and (ii) experiments typically provide 
spectra that are directly related to the response func- 
tions - e.g. the optical absorption cross section of a gas 
is proportional to the imaginary part of the dipole-dipole 
molecular polarizability. In contrast, in the strong-field 
regime, where perturbative treatments become cumber- 
some, one normally computes the particular response of 
the system to the perturbation of interest by directly 
propagating the TDKS equations in real time. 

The vast majority of TDDFT applications have ad- 
dressed the first-order response of the ground-state sys- 
tem to weak electric fields - which can provide the ab- 
sorption spectrum, the optically-allowed excitation ener- 
gies and oscillator strengths, etc. This can be performed 
by linearizing the TDKS equations in the frequency do- 
main and casting the result into matrix- eigenvalue form, 
or by propagating the same equations in real time ap- 
plying a sufficiently weak dipole perturbation. In any 
case, the response function computed in this manner will 



be that of the ground state. If we want to analyze a 
TAS experiment, the objective is to obtain the response 
of the excited states that are visited by the system as it 
is driven by the pump pulse (i.e. the response function 
of a system out of equilibrium). This extension will be 
treated in Section [HAj 

Likewise, TDDFT can be used to compute strong 
field non-linear photo-electron spectra of atoms and 
molecules, for example with the method recently devel- 
oped by some of us™ These spectra, however, are also 
characteristic of the ground state, although, as it will 



be shown in Section II B the methodology can be easily 
extended to tackle the pump-probe case (time-resolved 
photoelectron spectroscopy^) . 



A. Attosecond transient absorption spectroscopy 

When an electromagnetic pulse passes through a gas 
sample, the molecules polarize, and this polarization 
modifies the otherwise free propagation of light - one 
of the consequences being the partial absorption of it. 
In a dilute gas, assuming the electric dipole approxima- 
tion and a sufficiently weak pulse, the dipole-dipole lin- 
ear dynamical polarizability entirely determines the po- 
larization of the medium, and therefore the amount of 
absorption. This is usually understood at equilibrium: 
the gas is formed by molecules at thermal equilibrium 
(perhaps at sufficiently low temperature so that they all 
can be considered to be at their ground state), and the 
only light pulse present is that whose absorption we want 
to measure. 

In the pump-probe situation discussed here, however, 
one wishes to compute the absorption of a probe pulse 
by a set of molecules that is also irradiated by a pump, 
either simultaneously or with a given delay. The task is 
therefore to compute the response of the electric dipole 
with and without the probe pulse - the difference being 
the excess of polarization, responsible for the absorption 
of the probe. We will assume, as it is often the case, that 
the pump pulse is intense, whereas the probe is weak and 
can be treated in first order perturbation theory. 

This situation is amenable to a generalized definition of 
response functions, such as the one given in the appendix 
of Ref . [41] and also discussed in detail in Ref . |42j We will 
review here this definition, adapting it to the pump-probe 
situation. Let us depart from a Hamiltonian in the form 
(we only treat the electric part neglecting the magnetic 
term of the electromagnetic field): 



Ho[S](t)=H + S(t)V 



(5) 



where T~L is the static Hamiltonian that defines the sys- 
tem itself, and S(t)V is the coupling to the "pump" laser 
pulse. This is the "unperturbed" Hamiltonian, that con- 
tains only the pump pulse; the full Hamiltonian results 
of the addition of the probe pulse f(t)V: 



H(t)=H + S(t)V + f(t)V. 



(6) 
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The evolution of the system is given by: 



(7) 



and initially (t = to, some time before the arrival of both 
pump or probe), the system is at equilibrium: 



^5 Peq 



= 0. 



(8) 



For a fixed pump £, we may assume the system evolution 
to be a functional of the probe shape: p = p[f], and we 
may expand p in a Taylor series (in the functional sense) 
around / = 0: 



(9) 



n=0 



where po is the unperturbed system evolution (i.e. the 
evolution of the system in the presence of the pump 
only: / = 0), and p n is n-th order in the perturbation: 
p[Xf] = \ n p[f). The system response to this perturba- 
tion is measured in terms of the expectation value of an 
observable A, which can likewise be expanded: 



CO 

A(t) = Tr [p(t) A] =J2Mt)> 



n=0 



where 



A n (t) = Tr Pn (t)A 



(10) 



(11) 



For sufficiently weak probes, we are only interested in the 
first term: 



SA(t) = A(t) - A (t) « = Tr pi(t)i 



(12) 



which is linearly related to / through a pump-dependent 
response function: 

/CO 
dt' f(t')x A ^[S}(t,t') . (13) 
-co 

The response function is given by: 41 

XAyWJ) = ~ T± {pe, [Mm), V H [£](t')] } 

(14) 

Inside the commutator, the operators appear in the 
Heisenberg representation: 



d H [£}(t) = U[£}(to,t)dU[£}(t,t ), 



(15) 



where U[£] is the time propagation operator in the pres- 
ence of the pump only - hence the functional dependence 
on £ . We keep this functional dependence explicit in the 
notation for xa yl^KM')? to stress that it is a property 
of both the system (defined by the static Hamiltonian 



H) and of the pump shape £ , as opposed to the con- 
ventional response functions, which are only system de- 
pendent. Note also its dependence on two times t and 
t' ', which cannot be reduced to only one by making use 
of time-translational invariance, as it is customary when 
working at equilibrium. 

The response itself, 5A(t), will be a functional of both 
pump and probe pulses. If we take its Fourier transform 
we may write it as: 



5A[£J](uj) 



dt'f{t') XA v[£]{u,t'). (16) 



In order to compute the response function, one can use 
as a probe a delta perturbation, i.e. /(£') = \5(t' — r), 
which permits to identify: 

Xa #[S\{u,t) = ±5A[S,X5 t ]{u). (17) 

The action of such a delta-perturbation applied at the 
instant r on the system is given by: 



|$(t-> T +)> =e- iW |$(r)). 



(18) 



From now on we will restrict the discussion to pure states, 
since ensembles are not needed for the results that will 
be shown below. However it can easily be extended to 
general ensembles. We also restrict the discussion to a 
specific response function: the dipole-dipole polarizabil- 
ity a[£](M') ~ Xf) dI^K^'O' wnere both A and V are 
the atomic or molecular dipole operator D- taking into 
account that, for the frequencies that we are dealing with, 
the dipole of interest is that of the electrons, and the 
clamped nuclei approximation can be used. Moreover, 
we choose to work with light polarized in the x direction, 
so that: 



N 



(19) 



where N is the number of electrons. The expectation 
value of this electronic dipole is an explicit functional of 
the time-dependent density, and so is its variation: 



5D[£J](o 



d 3 r 5n(r, uo)x , 



(20) 



where 5n(r,u) is the Fourier-transformed difference be- 
tween the electronic densities obtained with and without 
the probe pulse. This straightforward formula in terms 
of the density is what makes TDDFT specially suited 
for these computations: we may safely utilize the Kohn- 
Sham system of non-interacting electrons. The delta per- 
turbation, Eq. (18), must be applied to each one of the 



Kohn-Sham orbitals, and takes the following form: 



<p(r,t- 



AXx 



¥i{r,r) . 



(21) 



The absorption of a particular probe pulse /, is deter- 
mined by the induced polarization, given by 5D[£, f](w). 



We will compute the dynamical polarizability a[£](cj,r), 
which is the polarization induced by a delta perturbation, 



<x[£](u>,t) = j5D[£,\5 t ](u>), (22) 

since it allows to compute any particular response 
through the integration of Eq. (16), as long as the probe 



is weak. In particular, we will look at the imaginary part 
of a[£](uj, r), which is the part responsible for absorption. 
Note, finally, that in a 3D situation the polarizability is 
not a scalar but a tensor, since there are three possible 
light polarization directions, and three components for 
the system dipole moment. In most cases, one is inter- 
ested in the trace of this tensor, an averaged quantity 
that corresponds to the absorption of a randomly ori- 
ented sample of molecules. 



B. Time- resolved photoelectron spectroscopy 




The photoelectron spectra presented in this work are 
produced within TDDFT using the recently introduced 
Mask Method. 40 This method is based on a geometrical 
partitioning and a mixed real- and momentum-space time 
evolution schemed In the following, we summarize the 
main traits of the technique (we refer the reader to Ref.l40l 
for a complete description), and demonstrate how it can 
be straightforwardly applied to the non-equilibrium sit- 
uation required by pump-probe experiments. 

In photoemission processes a light source focused on 
a sample transfers energy to the system. Depending on 
the light intensity electrons can absorb one or more pho- 
tons and escape from the sample due to the photoelectric 
effect. In experiments, electrons are detected and their 
momentum is measured. By repeating measurements on 
similarly prepared samples it is possible to estimate the 
probability to measure an electron with a given momen- 
tum. From a computational point of view, the descrip- 
tion of such processes for complex systems is a challeng- 
ing problem. The main difficulty arises from the necessity 
of describing properly electrons in the continuum. 

In typical experimental setups, detectors are situated 
far away from the sample and electrons overcoming the 
ionization barrier travel a long way before being detected. 
The distances, that electrons travel are usually orders 
of magnitude larger than the typical interaction length 
scales in the sample. During their journey towards the 
detector, and far away from the parent system, they prac- 
tically evolve as free particles driven by an external field. 
The solution of the Schrodinger equation for free elec- 
trons in a time dependent external field is known analyt- 
ically in terms of plane waves as Volkov states. It seems 
therefore a waste of resources to solve the Schrodinger 
equation numerically in the whole space if a considerable 
part of the wave function can be described analytically. 

In order to take advantage of the previous observations 
we partition the space according to the scheme in Fig. [I] 



FIG. 1. Schematic description of the space partition imple- 
mented by the mask method. A mask function (a) is used 
to implement the spatial partitions (b). In region A (inter- 
action region) the TDKS equations are numerically solved in 
real space while in B (free propagation region) electrons are 
evolved analytically as free particles in momentum space. Re- 
gion C is where if a and ips overlap. 

(b). The space is divided into two regions, A and B; 
the inner region A, containing the system with enough 
empty space around, is where electrons are allowed to 
interact with each other and with the system, and re- 
gion B, defined as the complement of A, is where elec- 
trons are non-interacting and freely propagating. Every 
KS orbital <Pi(r) can be decomposed accordingly with 
<fi(r) = tpf(r) + (pf(r), so that ipf(r) resides mainly in 
region A and <pf(r) mainly in region B. 

The geometrical partition is implemented by a smooth 
mask function M(r) defined to be one deep in the interior 
of A and zero outside (see Fig. [I] (a)): 

^(r) = M(r)^(r) , (23) 
ipfir) = (1 - M(r))^(r) . (24) 

The mask function takes care of the boundary conditions 
in A by forcing every function to be zero at the border. 
In order to give a good description of functions extending 
over the whole space it is convenient to represent ipf(v) 
orbitals in momentum space (pf{p). 

A mixed real and momentum-space time evolution 
scheme can then be easily derived following the geometri- 
cal splitting. Given a set of orbitals at time t their value 
at a successive time t + At is provided by 

f (^f(r, t + At) = M(r)e- i&At ipf (r, t) 

\ 0f (p, t + At) = e-'^^ 1 A Vf (P, t) + #(p, t + At) 

(25) 
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with H being the effective single-particle TDDFT Hamil- 
tonian, A(t) the total external time dependent vector 
potential (the coupling with the external field is conve- 
niently expressed in the velocity gauge), and 



<p?(p,t+At) 



(2tt) 3 / 2 



dr(l-M(r))e- lHA Vf (r, t)e 1 ^ 

(26) 

constituting the portion of electrons leaving the system 
at time t + At. At each iteration in the evolution the out- 
going components of y>f(r) are suppressed in the interac- 
tion region by the multiplication with M(r) while being 
collected as plane waves in <pf(p) via (pf(p). The re- 
sulting momentum space wavefunctions are then evolved 
analytically simply by a phase multiplication. 

The advantage of using such an approach resides in the 
fact that we can conveniently store the wavefunctions 
on a spatial grid inside A while treating wavefunctions 
in B (and therefore the tails extending to infinity) as 
free-electrons in momentum space. Moreover the mask 
function introduces a region C, where the wavefunctions 
in A and B overlap (see Fig. [I]) and that acts as matching 
layer. In spite of the fact that, from a theoretical point of 
view, the matching between inner and outer region could 
be performed on a single surface, from a numerical point 
of view, having a whole region to perform the matching 
is more stable and less influenced by different choices of 
spatial grids. 

From the momentum components of the orbitals in B it 
is possible to evaluate the momentum-resolved photoelec- 
tron probability distribution as a sum over the occupied 
orbitals 



P(p)« lim £>f(p,t)| s 



(27) 



the limit t — >• oo ensuring that all the ionized components 
are collected. This scheme is entirely non-perturbative; 
in a pump-probe setup, it does not assume linearity in 
either pump or probe. Therefore, it can be applied in 
the same manner when two pulses are present than with 
one pulse only, as it was shown in Ref. 40. Like in 
Sec. |II A| we can generalize the previous derivation to ad- 
dress transient photoelectron spectroscopy (spin-, angle- 
and energy-resolved) in practice by employing a pump- 
probe scheme and performing numerical simulations with 
two time delayed external pulses. A TRPES map is then 
generated by performing a computation for each different 
time delay. 

From P(p) several relevant quantities can be calcu- 
lated. The energy-resolved photoelectron probability 
P{E), usually referred to as photoelectron spectrum 
(PES), can be obtained by integrating P(p) over solid 
angles 



P(E = p 2 /2) 



<m p p(p) . 



(28) 



47T 



The angular- and energy-resolved photoelectron proba- 
bility P(#, </>, E), or photoelectron angular distribution 



(PAD), can easily be evaluated by expressing P(p) in 
polar coordinates with respect to a given azimuth axis. 
It is noteworthy that during the evolution defined in 



Eq. 25 the part of the density contained in A transferred 
to B is not allowed to return. Clearly, in cases where the 
external field is strong enough to produce electron orbits 
crossing the boundary of A and backscattering to the 
core the mask method provides a poor approximation. 
In these cases a bigger region A or a more refined scheme 
must be employed 4 ^. The laser fields employed in this 
work are weak enough that we can safely assume region 
A to always be sufficiently large to contain all the relevant 
electron trajectories. 



III. RESULTS 

The theory described in the previous section has been 
implemented in the octopus code^ 4 ; we refer the reader 
to Refs.l45land 46 for the essential points of the numerical 
methodology. 

In order to simplify the illustration of the results a 
clamped ion approximation has been used in the cal- 
culations for the molecular case. Further inclusion of 
the ionic motion could b e don e at the semi-classical level 
with Ehrenfest dynamic d 45 * 47 * - already implemented in 
the code - and without any essential modification to the 
theory presented. 



A. One-dimensional model Helium 

As first example we study the absorption spectrum of 
an excited one-dimensional soft-Coulomb Helium atom. 
This is an exactly solvable model that provides a useful 
benchmark to test different approximations. We will first 
discuss the exact solution, and later apply TDDFT. A 
more realistic 3D model will be presented in the next 
section. 

The ID model of the Helium atom is defined by the 
following Hamiltonian: 



H(t) = T + V ext (t) + V e( 



(29) 



where 

V ex t 



■£(t)( Xl +x 2 ) (30) 



is the external potential: the electron Coulomb interac- 
tion l/\x\ is softened to 1/Vl + x 2 . The coupling with 
the external time-dependent field £(t) is expressed in 
length gauge, and electrons are confined to move along 
the x direction only. Finally, 



2 [dxj + dxj 



(31) 
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FIG. 2. Out of equilibrium absorption spectrum as function 
of the pump laser frequency for one-dimensional Helium. The 
system is driven out of equilibrium by 45 cycle sin 2 envelope 
laser pulses of intensity I — 5.26 x 10 11 W/cm 2 , at differ- 
ent carrier frequencies and then probed right after. Maximal 
response is observed for frequencies close to the first optical 
transition uo-p = 0.533 a.u.. 

is the kinetic energy, and the electron-electron interaction 
is 



This model is numerically solvable given the exact map- 
ping discussed in Refs . l48l and [49| where it is proved that 
the many-body problem of N electrons in one dimension 
is equivalent to that of one electron in N dimensions. The 
wave functions and other necessary functions are repre- 
sented on a real space regular grid; a squared (linear for 
one-dimensional TDDFT) box of size L = 200 a.u. and 
spacing Ax = 0.2 a.u. has been employed in all the cal- 
culations. 

In order to illustrate how an external field can modify 
the optical properties of a system in Fig. [2] we show a scan 
of the non-equilibrium absorption spectrum generated by 
a 45 cycle sin 2 envelope pulse with intensity / = 5.26 x 
10 11 W/cm 2 at different angular frequencies 0.51 [a.u.] < 
uo-p < 0.59 [a.u.]. The plot displays Im[a[£ , t](cj, cj-p)], 
choosing r right at the end of the pump pulse. 

It can be seen how the absorption around the first ex- 
citation frequency 0.533 a.u. is strongly diminished when 
the frequency of the pump resonates with that frequency. 
In that situation, an absorption peak also appears around 
the excitation frequency corresponding to the transition 
from the first to the second excited state, in our case at 
0.076 a.u. This behavior is a direct consequence of the 
fact that the laser is pumping the system to the first ex- 
cited state and this process is more efficient for a field 
tuned with the excitation energy. The absorption spec- 
trum is therefore a mixture of the one corresponding to 
the unperturbed ground state and that of the first excited 
state. 
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FIG. 3. Comparison of absorption spectra calculated in dif- 
ferent approximations for a one-dimensional Helium model. 
The filled curves are the spectra for the unperturbed systems 
while the solid lines are the spectra of the system excited by a 
laser as in Fig. [2] resonant with the first allowed optical transi- 
tion: exact time-dependent Schrodinger equation uo-p — 0.533 
a.u. (in blue), EXX uj v = 0.549 a.u. (in red), and LDA 
(jj r = 0.475 a.u. (in green). The dashed blue line is the 
absorption of the system perturbed by a 180 cycle laser and 
probed at t = 30.62 fs, where the population on the excited 
state is maximal. The lines have been shifted by a vertical 
constant to facilitate the comparison between results. 



In order to analyze this point further, a cut at the res- 
onant frequency is displayed in the lower (blue) curves 
of Fig. [3j The filled curve represents the spectrum ob- 
tained from the system in its ground state while the solid 
line corresponds to the spectra of the system excited by 
a laser pulse with a frequency resonant with the first 
optical transition and probed after the perturbation. By 
direct comparison of the two spectra it is easy to discrim- 
inate the peaks associated with the ground-state absorp- 
tion from the ones characterizing the absorption from 
the excited state. In particular, the peaks related to 
the ground-state absorption are located at the energies 
co>o^i = 6i — eo = 0.533 a.u. corresponding to the transi- 
tion from the ground (eo = —2.238 a.u.) to the first ex- 
cited state (ei = —1.705) and uo^s = 63 — 60 = 0.672 a.u. 
corresponding to the third excited state (63 = —1.566 
a.u.) - the direct excitation of the second excited state 
is forbidden by symmetry. The solid curves show finger- 
prints of the population of the first excited state, namely 
the peaks corresponding to transitions from that first ex- 
cited state to others: in particular, the peak appearing at 
the low energy Ui^2 = £2 — £\ = 0.076 a.u. is associated 
with the transition from the first excited state c\ to the 
second one 62 = —1.629 a.u. 

These spectra contain information that is not con- 
tained in the equilibrium ones. For example, let us con- 
sider the spectra that would be produced by each sin- 
gle eigenstate, given by the state-dependent dynamical 
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polarizabilities, which may be written in the sum-over- 
states form as: 



00 — (Cj — €i) + i0+ 

w + (ej - ei) + i0+ 



(33) 



The poles of this function provide us with the eigenvalue 
differences Cj — e^; if this value is positive, the correspond- 
ing term is associated with a photon-absorption process; 
if it is negative, with a stimulated emission term. The 
weight associated with each one of these poles provide us 
with the dipole coupling matrix elements (^i\D\^j). 

During the time evolution, the wavefunction can be 
expanded on the basis of eigenstates of the unperturbed 
system Sfr(t) = rji(t)^i. When the system is probed at 
a certain time t, the resulting spectrum can be thought as 
a linear combination of the spectra produced by each sin- 
gle eigenstate. An analysis of the transient spectrum may 
therefore provide information about the mixing weights 
rji, and about excitation energies and dipole couplings 
between excited states - information that is absent in 
equilibrium ground-state linear response. 

In our case, we find by direct projection of the 
time-dependent wave function onto the eigenstates, that 
the system after the pulse is composed mainly of the 
ground and the first excited state with weights, 1 770 1 2 = 
|(*o|*(*)>| 2 = 0.7120 and |t?i| 2 = |(*i|tf(t))| 2 = 0.2876. 
The same information can be recovered by comparing 
the perturbed and the unperturbed spectrum at Eq-^i . 
At this energy we only have the contribution coming 
from ^0 ~- ^ ^1 and its inverse ^1 — >• ^o- The peak 
height of the perturbed spectrum after the laser pulse 
h t is therefore a combination of the heights associated 
with the ground ho and the excited hi states: h t = 
\Vo\ 2 ho + |^i| 2 ^i- At this energy ho = —hi due to the 
nature of the transition 1—^0 the ratio a — h t /ho = 
I ?7o 1 2 — 1 2 = 0.4258 thus gives direct information 
on the difference of the mixing weights. Complement- 
ing this information with a two-level system assumption 
\Vo\ 2 + \m\ 2 = 1 we obtain \f] \ 2 = (1 + a)/2 = 0.7129 
and |?7i| 2 = (1 — a)/ 2 = 0.2871 in good agreement with 
the results calculated by direct projection of the wave- 
function. 

In Fig. [4] we display the population weights for two 
different laser pulses. The red lines correspond to the 
same laser pulse as in Fig. [3] while the blue lines pertain 
to a four times longer laser with the same parameters 
(intensity, envelope shape and carrier frequency) and 180 
optical cycles. For both lasers the populations of both the 
ground and first excited state at each time almost sum 
to one, indicating an essential two- level dynamics. In the 
case of the long pulse we observe a maximum (minimum) 
of the population over the excited (ground) state at t = 
30.62 fs. This behavior can be understood in terms of 
Rabi oscillations. 




/ Is 



FIG. 4. Exact population on the ground |?7o(£)| 2 = 
|(^o|^(£))| 2 (solid lines) and first excited \r)i(t)\ 2 = 
|(^i|^(t))| 2 (dashed lines) states as a function of time for 
different laser pulses. In red a 45 cycles pulse with parame- 
ters as in Fig. [3] and in blue a longer 180 cycles pulse with 
the same parameters. 



A Rabi oscillation is a fluctuation behavior of states 
occupation occurring due to the interaction of an oscilla- 
tory optical field in resonance with a two-level system.^! 
The occupation probability alternates with the Rabi fre- 
quency Q(t) — f(t)/io->ii where fio^i is the dipole tran- 
sition matrix element between the states and f(t) is the 
electric field envelope. Extremal points of the popula- 
tions should be located at times where the pulse area 
®W = I-00 d T f(t) Vo-ti is an integer multiple of 7r, 
©(£) = U7T. With the numerically calculated matrix el- 
ement /jLo->i — 1-11 a.u. the first maximal population 
of the excited state is expected at t = 30.65 fs, in good 
agreement with what is observed. The absorption spec- 
trum at this time, as shown in Fig [3] (dashed blue line), 
displays a considerable enhancement at oj\^2 and a neg- 
ative emission peak at c^o^i as expected from a pure 
excited state. 

It is interesting to study the same model with TDDFT 
instead of with an exact treatment in order to address the 
performance of available (mainly static) xc-functionals. 
In Fig. [3] we display results obtained with TDDFT, em- 
ploying two different exchange-correlation (xc) functional 
approximations: exact exchange^ (EXX) in red, and 
one-dimensional local density approximation^ (LDA) in 
green. The calculations were performed in the adiabatic 
approximation using the same parameters as in the ex- 
act case. The laser frequency was tuned to match the 
first optical transition appearing at: uj-p = 0.549 a.u. for 
EXX, and uo v = 0.475 a.u. for LDA. 

The unperturbed spectrum (solid curve) provided by 
EXX is in good agreement with the exact calculation, 
and the perturbed one qualitatively reproduces the ex- 
act result. In particular the new peak appearing at low 
energy associated with the transition 1 — >• 2 is well repre- 
sented. In contrast LDA is only capable of reproducing 



FIG. 5. Comparison of the absorption spectra of unperturbed 
(filled curve) and perturbed He atom probed at r = 5.39 
(solid line) and after the end of the pulse r = 8.68 fs (dashed 
line). The spectrum range is below the ionization threshold. 
The atom is excited by a 45 cycle sin 2 envelope laser pulse 
polarized along the x-axis with carrier cj-p = 0.79 a.u. reso- 
nant with the first optical transition, intensity I = 2.6 x 10 12 
W/cm 2 . 



one peak both for the perturbed and unperturbed cases. 
This is due to the known problem of asymptotic exponen- 
tial decay of the functional that in this one-dimensional 
example supports only a single bound excited state. 

A common feature of both approximations is consti- 
tuted by the presence of negative values in the perturbed 
spectra. This can be tracked down to the lack of memory 
in the adiabatic xc-functional approximation 53 . The lack 
or wrong memory dependence in the functional results in 
slightly displaced absorption and emission peaks associ- 
ated with the same transition. This fact, analyzed in the 
light of Eq. [33j results, at the transition energy, in a sum 
of two Lorentzian curves with different sign and slightly 
different centers. This explains why we get two inverted 
peaks where we should have only a single one going from 
positive to negative strength as we populate the excited 
state - as shown by the exact (blue) curves in Fig. [3] 



B. Helium atom in 3D 

In this section we study the real Helium atom. We 
employed the EXX functional and discretized TDDFT 
equations on a spherical box of radius R = 14 a.u., spac- 
ing Ar = 0.4 a.u. and absorbing boundaries 2 a.u. wide. 

We begin by investigating the changes in absorption 
of He under the influence of an external UV laser field 
driving the system with the frequency of the first dipole- 
allowed excitation. To this end we used a 45 cycle sin 2 
laser pulse in velocity gauge with carrier uu-p = 0.79 a.u. 
resonant with the Is 2 —> ls2p transition, of intensity 
I = 2.6 x 10 12 W/cm 2 polarized along the x-axis. In 
Fig. [5] we show a comparison of the absorption spec- 



FIG. 6. Helium transient absorption spectrum scan for dif- 
ferent time delays. Pump laser pictured in the upper panel is 
the same as in Fig. [5] 



trum for the unperturbed atom (filled curve) and the per- 
turbed one probed with a delta perturbation right after 
the pump pulse at r = 8.68 fs (dashed line). The com- 
parison presents many traits similar to the ones discussed 
in Sec. MI Al for the one-dimensional Helium model. In 
particular, fingerprints of the population of the first ex- 
cited state can be observed in the appearance of a peak 
in the gap at uj2p^3s — 0.079 a.u. associated with the 
transition \s2p — » ls3s. The second peak, associated 
with the transition Is 2 — >• ls2p, presents height changes 
correlated with the former one. We also obtain the small 
artifacts, such as the energy shifts and the negative values 
attributed to the xc-kernel memory dependence discussed 
previously. 

Additional details on the excitation process can be ac- 
quired by expanding the time dimension of the absorp- 
tion spectrum. In Fig. [6] the time resolved absorption 
spectrum (TAS) map is displayed. The map was pro- 
duced by probing the system at different time delays. 
As the delay is increased we observe the build-up of the 
peak associated with the state being pumped by the laser 
pulse at uj2 P ^3s- This changes are reflected in the oscil- 
lations of the ground state first optical peak at ui s ^2 P - 
In TDDFT the knowledge of the wavefunction is lost in 
favor of the density, which does not allow us to do a pop- 
ulation analysis based on simple wave function projection 
The transient absorption spectrum, on the other hand, is 
an explicit density functional, and its computation with 
TDDFT may help us to understand the evolution of the 
state populations. 

The peak appearing in the gap presents a maximum 
at r = 5.39 fs that emerges before the end of the pump 
pulse (r = 8.68 fs). This peak is associated only with 
the transition from the ls2p — >• ls3s and therefore its 
height is proportional to the 2p excited state population. 
The oscillation can then be interpreted in terms of Rabi 



physics as discussed in Sec. Ill A 
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FIG. 7. Helium transient photoelectron spectrum in logarith- 
mic scale. The pump laser (upper panel) is the same as in 
Fig. [5] and the probe is a 40 cycles trapezoidal laser pulse 
with 8 cycles ramp, uo p = 1.8 a.u., I = 5.4 x 10 9 W/cm 2 
aligned with the pump pulse. 



Further insight can be achieved investigating the pho- 
toemission properties of the system. In Fig. [7] we show 
the TRPES map, as calculated in a pump-probe set up. 
Photoelectrons are calculated with the technique out- 
lined in Sec. |IIB| The pump pulse is the same as the 
one employed for TAS. The probe is a 40 cycles trape- 
zoidal laser pulse (8 cycles ramp) with carrier frequency 
oo p = 1.8 a.u., intensity I = 5.4 x 10 9 W/cm 2 , polar- 
ized along the x-axis and is weak enough to discard non- 
linear effects. We performed a scan for different time 
delays, measuring each delay as the difference from the 
probe center to the beginning of the pump. Negative 
delays correspond to the situation where the probe pre- 
cedes the pump. Moreover, in order to include all the 
relevant trajectories a spherical box of R = 30 a.u. was 
employed, and photoelectrons were recorded only during 
the up-time of the probe pulse. 

The TRPES map in Fig [7] shows three main features 
at Ei = 0.66 .a.u, E2 = 0.88 a.u. and E3 = 1.67 a.u.. 
In our case the probe pulse is weak, and photoelectrons 
escaping the system undergo photoelectric-effect energy 
conservation. A bound electron can absorb a single pho- 
ton and escape from the atom with a maximum kinetic 



E/a.u 



0/deg 
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energy E 



/p, where uj p is the probe carrier fre- 



quency and Ip is the field-free ionization energy. The 
ionization potential can be evaluated in DFT as the neg- 
ative energy of the highest occupied KS orbital (HOMO) 
Ip = £2s = 0-92 a.u. Thus the peak appearing at E 2 
is energetically compatible with photoelectrons emitted 
from the 2s level: E2 = oj p — Ip. Consistently, this peak 
is the only one appearing at negative delays where the 
pulses do not overlap. Moreover, the peak strength is 
weakly varying with the delay while slightly shifting to- 
wards lower values around 3 fs in accordance with TAS 
findings. At about the same delay time the peak at E3 



FIG. 8. Energy- and angular- resolved photoelectron spec- 
tra for Helium at fixed delay r = 8.99 fs. Panel (c) dis- 
plays a logarithmic scale PES P(E) comparison at fixed de- 
lays r = —1.69 fs (red) and r = 8.99 fs (green). The other 
panels depict normalized PADs P(0, 0, E) with polar coordi- 
nates referred to axis z at fixed delay r = 8.99 fs and energy: 
(a) Ex = 0.66 a.u., (b) E 2 = 0.88 a.u., and (d) E 3 = 1.67 a.u.. 
White crosses mark the intersection between the probe polar- 
ization axis and the cutting sphere. 



begins to emerge. This peak corresponds to emission 
from the pump-excited 2p state E3 = uj-p + uj p — Ip. 
It is a process where the atom, initially in the ground 
state, absorbs a photon from the pump and gets excited 
to the 2p bound state. The subsequent absorption of a 
probe photon frees the electron into the continuum. The 
peak at E\ is understood in terms of pump photons only: 
E3 = 2cup — Ip. The ionization mechanism shares the 
first step with the E3 process, namely a 2s — » 2p excita- 
tion produced by the absorption of a uop photon. In the 
second step the electron is directly excited to a contin- 
uum state by the absorption of a second up photon. In 
the linear regime, the direct photoionization cross-section 
decays exponentially with energ 5P. For this reason and 
due to the disparity in intensity between pump and probe 
this ionization channel is by far the most favorable one. 

In direct photoemission processes, the photoelectron 
angular distribution (PAD) contains information about 
the electronic configuration of the ionized state. 55 In or- 
der to support the energetic arguments PADS P(0, (/>, E) 
at r = 8.99 fs are presented in Fig [8] (a), (b), (d) together 
with cuts on TRPES map at r = fs and r = 8.99 fs (b). 
For each energy marked in Fig. [8] (c) we perform spheri- 
cal cuts of the photoemission probability on energy shells 
at E — E\, £?2, E3. Each cut is then plotted in polar co- 
ordinates with being the angle from the z-axis and <p 
the angle in the x?/-plane measured from the x-axis. In- 
tersection of the lasers polarization axis with the sphere 
are marked with a white cross. 

Fig [8] shows clearly that photoelectrons at E\ (a) and 
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E% (d) have similar nature compared to E 2 (b), in agree- 
ment with the energy analysis. Electrons emerging with 
a kinetic energy of E2 are emitted from a 2s state, and 
symmetry of the orbital is imprinted in the photoelec- 
trons angular distribution. In order to understand the 
PAD features it must be taken into account that 2s elec- 
trons are perturbed by a laser with a specific polarization 
direction that breaks the rotational symmetry. The laser 
transfers maximal kinetic energy along directions parallel 
with the polarization and minimal along the perpendic- 
ular plane and, if non-liner effects can be discarded, it 
induces a geometrical factor of the form | A • p| where A 
is the polarization direction and p the electron momen- 
tum. For this reason electron emitted along (j) = 90° , and 
270° are strongly suppressed, and panel (b) is compatible 
with the spherical symmetry of a 2s state. 

In panel (a) electrons are excited to a p state and 
then ejected into the continuum by the absorption of 
two pump photons. The PAD displays marked emission 
maxima for the direction aligned with the laser polar- 
ization (indicated by white crosses). The extension in 
is narrower compared with the 2s emission in panel 
(b) consistently with ionization from a p x orbital. Of 
the three degenerate p orbitals the p x is the one pro- 
ducing the strongest response. Signatures of p y and p z 
response can be identified in the non vanishing PAD on 
the 7/z-plane around (j) = 90°, 270°. Such perpendicu- 
lar response indicates a degree of non-linearity induced 
by the pump. Similar considerations hold for panel (b) 
where the p state excited by the pump is probed with u p . 
As before the emission is mainly from a p x state. 



C. Ethylene molecule 

In this section we extend our calculations to the treat- 
ment of the Ethylene molecule (C2H4) and show how 
these techniques permit to study the time-dependence of 
molecular electronic states. In particular we report on 
the clear observation of a strong ir — » 7r* transition. 

In order to have a good description of states 
close to the ionization threshold we employed the 
asymptotically correct LB94 xc-functional in the adi- 
abatic approximation^. We choose the molecular 
plane to be in the x?/-plane with carbon atoms at 
coordinates (±1.26517,0,0) a.u. and hydrogens at 
(±2.33230,1.75518,0) a.u., (±2.33230,-1.75518,0) a.u.. 
The ion positions are held fixed during the time evo- 
lution. Norm-conserving Trouiller-Martin pesudopoten- 
tials are employed to describe core electrons of Carbon. 
Moreover TDDFT equations are numerically integrated 
on a spherical grid with spacing Ar = 0.3 a.u., radius 
R = 16 a.u. and 2 a.u. wide absorbing boundaries. 

We perturb the system with a 15 cycles (3 cycles 
ramp) trapezoidal laser pulse having carrier frequency 
uj v = 0.297 a.u. and intensity I = 1.38 x 10 11 W/cm 2 
polarized in the x-axis. The laser frequency and the po- 
larization direction are suited to excite mainly the molec- 
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FIG. 9. Comparison of the absorption spectra of unperturbed 
(filled curve) and perturbed Ethylene molecule (solid line) 
below the ionization threshold. The molecule is excited by a 
45 cycle sin 2 envelope laser pulse polarized along the x-axis 



with carrier uo-p 
W/cm 2 . 



0.297 a.u. of intensity I = 1.38 X 10 1 




FIG. 10. Ethylene molecule TAS. Pump laser pictured in the 
upper panel is the same as in Fig. [5] 



ular 7r — ^ 7r* transition. The absorption spectrum of the 
excited molecule probed after the pulse Fig. [9] shows the 
emergence of a peak associated with a population of the 
7r* state. Optical transitions from this excited state to 
high lying bound states occur at energies lower than the 
HOMO-LUMO gap as illustrated by the scheme in Fig. [9) 
Effects of the lack of memory in our xc-potential can be 
observed in the shifts of the peaks with respect to the 
known excitations of the unperturbed system. The char- 
acteristic excitations of a many-body system should not 
depend on the perturbation, unless we are in a strong 
light-matter coupling regime. 

The build up of the transient spectrum as a function of 
time is shown in Fig. 



He discussed in Sec I 
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In comparison with the case of 
IB the TAS map does not display 



any maxima during the pump time lapse due to the enve- 
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FIG. 11. Logarithmic scale TRPES for C 2 H 4 . The molecule 
is probed at different delays with is a 40 cycles trapezoidal 
laser pulse with 8 cycles ramp, oj p = 1.8 a.u., / = 5.4 x 
10 9 W/cm 2 polarized along the z-axis perpendicular with the 



pump. Pump laser (upper panel) is the same as in Fig. 10 



lope area not crossing tt by the end of the pulse. A pulse 
with larger area would reveal the firs Rabi oscillation. 

The TRPES map is presented in Figjll] Calculations 
were performed in a box of radius R = 30 a.u. and the 
probe pulse i s equal to the one used for the Helium atom 
in Sec. IIIB 1 namely a 40 cycles (8 cycles ramp) trape- 
zoidal pump at oo p = 1.8 a.u. and / = 5.4 x 10 9 W/cm 2 , 
but polarized along the z-axis. The choice of the polar- 
ization direction is important since, as we shall show, the 
spectra may reveal geometrical features that depend on 
it. 

The TRPES displays a behavior similar to the one dis- 
cussed in Sec. IIIIBI A set of constant crests can be ob- 
served in the energy range from 0.8 a.u. to 1.5 a.u.. These 
are the peaks associated with electrons residing in the 
ground- state and ejected by the probe pulse. In partic- 
ular the peak at E 2 = 1.37 a.u. corresponds to emission 
from the tt HOMO E 2 = oj p — e^, with = 0.43 a.u.. 
The 7r orbital is localized on the Carbon atoms with two 
density lobes lying in the xz-plane and nodes in the xy- 
plane (refer to Fig. 12 (a) for geometrical visualization). 
A probe laser orientation along y should suppress elec- 
trons on xy-pldine perpendicular to its polarization and 
therefore PAD P(9,^,E 2 ) in Fig. [12] (c) for = 0° is 
diminished also due to geometrical reasons. Signature of 
a 7r symmetry can be clearly observed in the oscillations 
with cj) presenting maxima at <fi = 90°, and 270° along the 
plane perpendicular to x that indicates a concentration 
along the C — C bond axis, and minima for = 0°, 360° 
consistent with a depletion in direction of each carbon 
atom. 

Separated by a probe photon uj<p at E% = uo p + uo-p — 
= 1.67 a.u. we find photoelectrons ejected from the 
7r* state. The intensity of the peak steadily increases in 
time accordingly to what is observed with TAS. Com- 
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FIG. 12. Angular- and energy- resolved photoelectron spectra 
for C2H4 at two fixed delay times. Panel (a) geometry of the 
process: p indicates photoelectron direction, A is the pump 
polarization vector, and A p the probe one. Panel (d), loga- 
rithmic scale PES P(E) for r = -1.69 fs (red) and r = 3.63 fs 
(green). Other panels depict normalized PADs P(9,<fi,E) at 
r = 3.63 fs and energies marked in (d): (b) E\ — 0.16 a.u., (c) 
Ei — 1.37 a.u., and (e) E3 = 1.67 a.u.. White marks indicate 
the position of the probe polarization vector (c), (e) (at the 
corners) and the pump one (b) on the sphere. 



pared with 7r, the 7r* orbital presents additional nodes on 
the plane perpendicular to the molecular bond and a field 
polarized along z is sensible to this kind of geometry. The 
PAD P(0, 0, £3) in Fig. 12 (e) displays strong suppres- 
sion of electrons along the ?/z-plane at <p = 90°, 270° and 
therefore presents a clear manifestation of photoemission 
from a 7r* state. 

Slow electrons ejected at E\ — 2ujj> — = 0.16 a.u. 
gradually increase and become the predominant ioniza- 
tion channel. The emergence in time of multi-photon 
peaks separated by oj-p indicates that the pump is strong 
enough to trigger non-linear effects. These electrons 
are ejected after the simultaneous absorption of pump 
photons. Electrons at E\ reach the continuum with an 
uj-p photon after the molecule has been excited to a 7r* 
state by another ujp photon. PAD should therefore carry 
again signs of 7r* symmetry. It must be noted, that in 
this case, 7r* electrons are probed with the pump itself, 
and therefore the laser polarization is along x. As al- 
ready discussed in the previous section, the laser polar- 
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ization carries a geometrical factor of the form A • p with 
A = Ax, that introduces a suppression along the yz- 
plane (</> = 90°, 270°). Unfortunately, this plane is pre- 
cisely where the 7r* photoemission minima should lay. For 
this reason PAD P(0,</),E 3 ) in Fig. [12] (b) is not suited 
to clearly discern a tt from a 7r* symmetry, and the sup- 
pression for = 0° along the x?/-plane is compatible with 
both structures. 



IV. CONCLUSIONS 

In this work we studied the problem of describing ul- 
trafast (attosecond scale) time-resolved absorption and 
photoemission in finite systems with TDDFT. We pre- 
sented the theory and discussed how it can be imple- 
mented, such that TDDFT can be successfully employed 
in the task of describing the dynamics of electronic ex- 
cited states in atoms and molecules. We illustrated the 
theory with three applications: one-dimensional Helium 
model, three-dimensional Helium atom, and Ethylene 
molecule. 

We studied the one-dimensional Helium atom per- 
turbed by an external time dependent field exactly, by 
solving TDSE. We showed how it is possible to recover 
information about state populations through a compari- 
son of the perturbed and unperturbed absorption cross- 
sections, and that the population evolution in time can 
be described in terms of Rabi physics. We then per- 
formed TDDFT calculations on the same model, and we 
may conclude that the results obtained with the EXX 
potential are in good agreement with the exact solution, 
although small artifacts appear due to the wrong descrip- 
tion of the functional memory dependence. 

Furthermore we investigated the Helium atom in a 
more realistic three-dimensional treatment using the 
EXX functional. We performed resonant pump-probe 
calculations monitoring both absorption and photoemis- 
sion properties of the excited atom. TAS turned out to 
be a sensible tool to monitor the build-up of the excited 
state, allowing to observe Rabi oscillations as a function 
of the time delay between pump and probe. TRPES also 
allowed the characterization of the excitation process in 
time. However, due to a dominant ionization channel 
associated with sequential two (pump) photons absorp- 
tion, the information about the excited state population 
was less apparent. Nonetheless PAD, being an observable 
sensitive to the geometrical arrangement of the ionized 
state, is a useful tool to discern the nature of each photo- 
electron peak. As a final example we considered the case 
of the Ethylene molecule, to study the time evolution of 
a 7r — >• 7r* transition. PAD for ejected electrons offered 
clear evidence that the states taking part in the process 
were indeed of tt and 7r* nature. 

The theoretical framework that we have developed is 
a useful tool to understand and control non-equilibrium 
electronic dynamical processes in nanostructures and ex- 
tended systems. New emergent properties of matter in 



the strong-coupling regime could appear that might give 
rise to new technological developments. Furthermore, 
monitoring electron and ion-dynamics provides funda- 
mental insights into structure (i.e. time-resolved crys- 
tallography) and chemical processes in biology and ma- 
terials science (e.g. for energy applications). There is 
plenty of room for new and fascinating discoveries about 
the behavior of matter under out-of-equilibrium condi- 
tions. 

Still, from the fundamental point of view, there is a 
clear need for the development of non-adibatic exchange 
and correlation functionals able to provide a reliable de- 
scription of non-equilibrium processes and strong-light- 
matter interaction. Clearly, the methods presented here 
will automatically benefit from any theoretical advance 
in this direction. Conversely, the developers of new func- 
tionals may take into account the correct description of 
pump-probe experiments as a useful quality test. 
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